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Abstract 

There are two promising scenarios that explain the ankle, which is a dip in the 
spectrum of cosmic rays at ~ 10 19 eV. A scenario interprets the ankle as the transi- 
tion from Galactic to extragalactic cosmic rays {ankle-transition scenario), while the 
other is that the dip caused by pair production on the cosmic microwave background 
radiation {proton-dip scenario). In this paper, we consider whether cosmogenic neu- 
trinos can be a clue to judge which scenario is favored. We calculated the fluxes 
of cosmogenic neutrinos following these scenarios with plausible physical parameter 
sets, and found several important features as follows. First of all, the neutrino flux 
at ~ 10 20 eV becomes much higher in the ankle-transition scenario as long as the 
maximum energy of the cosmic rays at sources is sufficiently high. On the other 
hand, the neutrino spectrum has a characteristic peak at ~ 10 16 eV in the proton- 
dip scenario on the condition that extragalactic protons significantly contribute to 
the observed cosmic rays down to 10 17 eV. Thus, we conclude cosmogenic neutrinos 
should give us a clue to judge which scenario is favored, unless these features are 
masked by the neutrino background coming from possible, powerful neutrino sources 
such as AGNs and GRBs. We also found an interesting feature that the neutrino 
flux at ~ 10 18 eV depends only on the cosmological evolution of the cosmic ray 
sources. That means cosmogenic neutrinos with the energy bring us information on 
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the cosmological evolution of the sources of ultra- high energy cosmic rays. Finally, 
we compare the fluxes of cosmogenic neutrinos with the expected sensitivity curves 
of several neutrino detectors, and conclude the detection of cosmogenic neutrinos 
in the near future is promising. 
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1 Introduction 

What energy is the transition point of Galactic cosmic rays (GCRs) and extra- 
galactic cosmic rays (EGCRs) in cosmic ray spectrum is an intriguing problem 
in cosmic ray astrophysics. The observed spectrum of cosmic rays, over more 
than 12 orders of magnitude in energy, can be described by a power-law shape 
with several spectral breaks [1]. The spectral breaks are interpreted as the 
transition points of sources to mainly contribute to the observed flux of cos- 
mic rays. A spectral break at ~ 10 15 ' 5 eV, so-called knee, has been thought of 
as the appearance of a maximum acceleration energy of protons at energetic 
objects in our Galaxy, like supernova remnants. 

At ~ 10 19 eV, a spectral dip, called ankle, has been observed. Traditionally, 
the ankle has been interpreted as a transition point from GCRs with a steep 
spectrum (oc E~ 3A ) to EGCRs with a harder spectrum (oc E~ 2 - 2 - 3 ). This is 
partly because cosmic rays above 10 19 eV could not be confined in the Galaxy 
by the Galactic magnetic field and their arrival distribution is highly isotropic. 
We call this traditional scenario ankle-transition scenario throughout this pa- 
per. In this scenario, it is an open problem how GCRs are accelerated up to 
10 19 eV in our Galaxy. On the other hand, it has been suggested recently that 
extragalactic protons with a steep injection spectrum (oc E~ 2 - 6 2 7 ) can repro- 
duce the ankle as a spectral dip due to Bethe-Heitler pair-creation with cosmic 
microwave background (CMB) photons [2,3]. This scenario, called proton-dip 
scenario throughout this paper, requires extragalactic protons at least down 
to 10 18 eV. In the proton-dip scenario, the second knee, which is another spec- 
tral break at ~ 10 17 7 eV, is interpreted as a transition point from GCRs to 
EGCRs. It is an intriguing problem for the determination of cosmic ray origin 
which scenario is favored. 

It is difficult to distinguish the two scenarios based on the energy spectrum 
of cosmic rays in observations. One of the key clues to determine the tran- 
sition point is cosmic ray composition. In the ankle-transition scenario, the 
composition of GCRs at ~ 10 19 eV is expected to be dominated by heavy nu- 
clei because of their capability of accelerating up to such high energy, while 
proton-dominated composition is required in the proton-dip scenario. Thus, 
composition measurements may provide us with useful information on the 
transition scenarios. However, the results of composition measurements are 
difficult to constrain transition scenarios because the uncertainty of hadronic 
interaction models in the simulations of extensive air showers obstructs the 
accurate determination of the composition [4,5]. 
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In this study, we focus on another key element, cosmogenic neutrinos, which 
are generated by photopion production between ultra-high-energy cosmic rays 
(UHECRs) propagating in intergalactic space and cosmic background pho- 
tons, and through the successive decay of producing pions and muons. Cos- 
mogenic neutrinos are a definite probe of EGCRs. The estimations of the flux 
of cosmogenic neutrinos have been performed since the discovery of the CMB 
[6,7,8,9,10]. A recent work of Ref. [11] pointed out that cosmogenic neutrinos, 
coupled with UHECR results, would provide a sufficient description of the 
properties of UHECR sources. 

Cosmic background radiation with the energy higher than the CMB (e.g., 
infrared (IR), optical, ultraviolet (UV)) can significantly contribute to the 
total flux of cosmogenic neutrinos because there is a large number of UHE- 
CRs which can generate neutrinos by interactions with such higher energy 
photons though the number of such photons is much smaller than the CMB 
[12,13,14]. Recent progress in the observations of high-redshift Universe al- 
lows constructing detailed models of the spectral energy distribution (SED) 
of IR to UV background (IR/UV below) radiation [15,16]. Refs. [13,14] dis- 
cussed the neutrino flux using the SED model constructed by Ref. [16]. Ref. 
[13] considered the propagation of UHE protons and estimated the resultant 
flux of cosmogenic neutrinos. Based on this work, it was suggested in Ref. 
[17] that cosmogenic neutrinos might be a key clue for determining the transi- 
tion energy, but detailed discussions were not performed. On the other hand, 
the authors of Ref. [14] estimated the fluxes of cosmogenic neutrinos on the 
assumptions of both a pure proton and mixed composition models. 

Motivated by these studies, we discuss whether cosmogenic neutrinos can be 
a clue to judge which transition scenario is favored in detail in this study. 
We calculate the spectra of cosmogenic neutrinos following the two transi- 
tion scenarios for this purpose for plausible physical parameter sets. Since the 
neutrino fluxes depend on several physical parameters, like the maximum ac- 
celeration energy of UHECRs, the minimum energy of EGCRs, the spectral 
shape of UHECRs, and the cosmological evolution of UHECR sources, we also 
investigate the parameter dependence of the neutrino fluxes and the capability 
of distinguishing between the two transition scenarios. We adopt an IR/UV 
background model other than a model used in Refs. [12,13,14]. The neutrino 
fluxes are normalized by the fluxes and spectral shapes of accompanying UHE- 
CRs obtained in observations. For simplicity, extragalactic magnetic field is 
neglected. 

The composition of EGCRs are assumed to be purely protons in this study. 
EGCR composition is an interesting problem in itself. The High Resolution 
Fly's Eye (HiRes) reported that the composition of cosmic rays above 10 18 eV 
is dominated by protons as a result of (V max ), the averaged depth of the 
shower maximum, measurement [4]. A recent result by Pierre Auger Observa- 
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tory (PAO) is consistent with the HiRes result within systematic uncertainty 
[5]. On the other hand, studies of muon content in the extensive air shower, 
another observable for UHECR composition, indicate a significant fraction of 
heavy nuclei above 10 19 eV [18,19]. The accurate interpretation of these com- 
position measurements is difficult because of our poor knowledge of hadronic 
interactions at ultra-high-energy, as mentioned above. The PAO also reports 
the positional correlation between the arrival directions of the highest en- 
ergy cosmic rays above 5.7 x 10 19 eV and nearby active galactic nuclei (AGNs) 
[20,21]. If these AGNs are really the sources of the observed events, this fact 
implies the highest energy cosmic rays are dominated by protons because of 
small deflections by the Galactic magnetic field. This implication is indepen- 
dent of the composition measurements. Thus, an assumption that EGCRs are 
dominated by protons up to the highest energies is reasonable. 

This paper is organized as follows: in Section 2, we explain our calculation 
method of cosmogenic neutrino fluxes in detail. In Section 3, we address our 
results and discuss the detectability of cosmogenic neutrinos taking neutrino 
oscillation into account. In Section 4, several uncertainties on the neutrino 
flux are discussed. We conclude in Section 5. 



2 Our Calculation Method 

In this section, our calculation method of the flux of cosmogenic neutrinos is 
explained. The neutrino flux can be obtained by calculating the number of 
neutrinos produced by propagating protons injected from a source, and then 
by integrating such neutrinos over all sources. We explain the propagation 
of UHE protons and their interactions with cosmic background photons in 
Section 2.1. In Section 2.2, our treatment of neutrino production is described. 
Then, we represent UHECR source models to calculate the total neutrino flux 
in Section 2.3. 



2. 1 Propagation of UHE Protons 

Protons propagating in intergalactic space interact with cosmic background 
photons and lose their energies through particle productions [9,22]. They also 
lose their energies adiabatically due to the cosmic expansion since they prop- 
agate over cosmological distance. We consider two interactions with cosmic 
background photons: photopion production and Bethe-Heitler pair creation, 
and adiabatic energy-loss. 

For cosmic background photons, we take into account not only the CMB 
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photons but also IR/UV background photons. Background photons with the 
energies higher than the CMB photons allow protons with the energies lower 
than the threshold of photopion production with the CMB to generate neu- 
trinos. Thus, they increase the neutrino flux significantly while they do not 
change UHECR spectrum because of the smaller number density of IR/UV 
photons compared with the CMB [13]. We adopt the best-fit model of the 
SED of IR/UV background photons calculated in Ref.[15]. 

Photopion production is the most essential process in this study. This process 
is treated stochastically. The mean free path of photopion production for a 
proton with energy E, in an isotropic photon field, \(E,z), is calculated as 
[23] 

1 1 °° 1 dn Smax 

e th mi n 



where s min = (ra^o + m p ) 2 c 4 , s max = m p 2 c A + 2Ee(l + (3), e th = (s min - 
m p 2 c 4 ) [2E(1 + /3)]" 1 . nip, m no , c, /3, e, and s are the proton mass, the neutral 
pion mass, the speed of light, the velocity of the proton in the unit of the speed 
of light, the energy of cosmic background photons, and the Lorentz invariant 
energy squared, respectively. a(s) is the total cross section of photopion pro- 
duction and dn^/de is the differential number density of cosmic background 
photons. a(s) is calculated by GEANT4, a Monte-Carlo simulation tool-kit 
which can simulate photomeson productions [24]. The GEANT4 can well re- 
produce the experimental total cross section of photopion production. 

The mean free paths in intergalactic space are calculated and used for every 
Az = 0.1 up to z = 5. The mean free paths at several redshifts z are shown 
in Fig.l. The mean free path at z = rapidly decreases above 8 x 10 19 eV 
where the channel of interactions with the CMB opens, which makes sharply 
cosmic ray spectrum steepening, so-called Greisen-Zatsepin-Kuz'min (GZK) 
steepening [25,26]. At z — 0, the minimum length of the mean free path A mm (,z) 
is about 4 Mpc. The interaction points of propagating protons are determined 
by a method used in Ref.[27], based on X(E,z) and A mm (2:). 

The inelasticity of photopion production, K(s), is approximated by a conven- 
tional method as K(s) = [1 — (itlcr 2 — m 7T 2 )/s] /2 for single pion production, 
where tocr and m w are the mass of the cosmic ray after an interaction (proton 
or neutron), the pion mass generated in the reaction, respectively. For mul- 
tipion production (see Section 2.2), we approximately adopt this formula by 
replacing m 7r with the total mass of pions. This assumption does not almost 
affect the flux of cosmogenic neutrinos shown in Section 3. A typical inelas- 
ticity is 0.23 at s ~ 1.6 GeV 2 , which is the Lorentz invariant energy squared 
at the delta resonance in the total cross section of photopion production. 
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Fig. 1. Mean free paths of protons for photopion production in cosmic background 
radiation including IR/UV photons (solid lines) and only CMB photons (dashed 
lines) at six redshifts. The energy-loss lengths of Bethe-Heitler pair creation of 
protons in cosmic background radiation from the CMB to IR/UV are also shown 
(dot-dashed lines). 



Cosmic background photons also make protons lose their energies through 
Bethe-Heitler pair creation, jry — > pe + e~. Fig.l also shows the energy- loss 
lengths of this process which are calculated following an analytical fitting func- 
tion given by Ref. [28]. This energy-loss process can be treated in a continuous 
energy-loss approximation since its inelasticity is very small (2m e /m p ~ 1CT 3 ). 
Since the energy-loss length is much shorter than the interaction length of 
photopion production for IR/UV photons at high-redshift Universe, neutrino 
production related to IR/UV photons mainly occurs at about less than z=l. 

The energy-loss rate of protons with energy E due to the cosmic expansion is 
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expressed as 



dE _ _a E _ 
dt a 



-H 



1/2 



E, 



(2) 



where ACDM cosmology with Vt m = 0.3, fi A = 0.7, and H = 71km s _1 Mpc -1 
is assumed. We also treat this process as continuous process during propaga- 
tion. 

The spectra of UHE protons and cosmogenic neutrinos are calculated as fol- 
lows: first, we consider protons with 10 16 -10 22 eV ejected from a source. This 
energy range is divided into 60 bins with A(\og 10 E) = 0.1 and 2000 protons 
are ejected in each bin initially regardless of UHECR injection spectrum. The 
propagation of protons is calculated taking the energy-loss processes and pion 
production into account, and then the resultant spectra of UHE protons and 
neutrinos at the Earth are recorded for each initial energy bin. The calcu- 
lation of propagation stops when the energies of propagating protons reach 
10 15 eV. Neutrinos are assumed to lose their energies only due to the adiabatic 
energy loss. Then, the resultant spectra sum up over all initial energy bins 
weighted by a given injection spectrum to obtain the spectra of UHE protons 
and cosmogenic neutrinos from a source. Finally, the total spectra of protons 
and cosmogenic neutrinos are calculated by integrating their spectra from a 
source over all sources taking UHECR source-evolution models into account. 
The injection spectra and UHECR source-evolution models adopted in this 
study are explained in Section 2.3. 



2.2 Neutrino Production 



Cosmogenic neutrinos are produced by the decay of charged pions generated 
by photopion production of propagating protons with cosmic background pho- 
tons. A charged pion decays into a muon and muon neutrino, ti + — > /U + z/ M or 
7r~ — > yTv^. The muon decays into an electron and neutrinos, /i~ — > e _ z/ e z/^ 
or fi + -> e + u e i7^ 

There are many reaction modes for pions in photopion production, since 
protons with sufficient energies could produce more than one pions. The 
GEANT4, which is used to calculate the total cross section, can also cal- 
culate exclusive cross section of each reaction mode, but it has some problems 
in parametrization. Therefore, we combine the experimental data of exclusive 
cross sections similarly to Ref.[29]. In addition, it takes a large CPU time to 
fully simulate photopion production one by one. Hence, we shall take a simpler 
treatment to calculate neutrino generation, as explained below. 

We consider pj — * pn , nn + as single pion production processes, and pj — > 
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£>7r + 7r~, n7r°7r + , pn°Ti as double pion production processes. We adopt exper- 
imental results of the cross sections of these reactions below s ~ 3 GeV 2 like 
Refs. [29,30,31]. These cross sections are not confirmed well by experiments 
at s > 3 GeV 2 . Therefore, we simply extrapolate the cross sections from 
s ~ 3 GeV 2 to higher energy. Then, we regard the total cross section minus 
the contributions of single and double pion production as the cross section of 
triple pion production, which generates 7r°, 7r + and 7r~. The above treatment 
is sufficient unless multi-pion production processes are significant. 

Once we judge that a propagating proton raises photopion production in cal- 
culation, the counterpart photon energy e, and s are determined by a method 
used in Ref.[23]. Next, it is determined which reaction is realized based on the 
probabilities proportional to the cross section of each reaction mode at s. The 
total energy of the produced pions is given by K(s)E. K(s)E is the energy 
of the pion for single pion production. For double or triple pion production, 
we assume the total energy, K(s)E, to be divided equally into all pions. Gen- 
erated charged pions decay into muons and muon neutrinos. The energies of 
the products are calculated using the two body decay algorithm. The ener- 
gies of the products of muon decay are calculated using the three body decay 
algorithm. The pions and muons decay immediately. 

Neutrinos are also produced by the beta decay of neutrons, n — » pe~v e , 
which result from charged pion production. These neutrinos can contribute 
to the observed neutrino flux if the mean free path of the beta decay, 7cr = 
0.92(E/10 20 eV) Mpc, is shorter than that of photopion production of a neu- 
tron which is assumed to be equal to that of a proton. At z — 0, for example, 
the mean free path of the photopion production is comparable with that of 
beta decay of a neutron with 10 20 6 eV. Hence, the neutron beta decay is also 
an important process for neutrino productions. This process generates only 
anti-electron neutrinos. The energies of the products are calculated using the 
three body decay algorithm. 



2.3 UHECR Source Models 



The total flux of cosmogenic neutrinos is calculated by integrating a neutrino 
flux from a source over all sources. The cosmological evolution of the number 
density and luminosity density of UHECR sources strongly affects on the 
flux of cosmogenic neutrinos. In this study, we consider four source-evolution 
models in which sources are uniformly distributed with cosmological evolution 
up to z = 5. 

The first is a uniform distribution without cosmological source-evolution (UNF) 
The second is a model following the luminosity density evolution of quasars 
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(QSO), which is parametrized by Ref. [32] 



/qso(^) oc 



(l + z) 3 (z<1.3) 

Constant (1.3 < z < 2.7) . (3) 
exp (l - ^) (2.7 < z) 



The third follows a star formation rate (SFR) deduced from the reconstruction 
of a IR/UV background flux used in this study, 



/sfr(z) oc < 



(l + 2; ) 3 - 5 z<1.2 

(4) 

;i + z)" L2 (1.2 <z) 



The fourth, the last model is motivated by an assumption that the observed 
UHECRs come from gamma-ray bursts (GRBs) [33]. Although the redshift 
distribution of GRBs is still controversial and it is different between pre-Swift 
and Swift bursts (e.g., [34]), one plausible rate history is that the GRB rate is 
more enhanced at higher redshifts due to metalicity effect [35]. It infers that the 
progenitors of GRBs would favor metal-poor stars, and host galaxies of GRBs 
would tend to have low metalicities [36,37]. In this paper, we shall adopt such 
an evolution model, and refer this GRB rate model as a " GRB-metalicity anti- 
correlation" model. More specifically, we shall assume a following rate history, 
/grb(z) oc (1 + z) 1A f SFR (z) , which is used in Ref. [38] based on the calculation 
of Ref. [39] . We use this parametrization in our fourth source-evolution model 
(SFR+GRBMAC below). As for the star formation rate, /sfr(z) in the third 
model is adopted. Our results do not depend on the local GRB rate which 
is somewhat uncertain so far, because we use the observed UHECR flux for 
the normalization of the flux of cosmogenic neutrinos. However, note that the 
lower local rate requires higher baryon loading per GRB if UHECRs come 
from GRBs. 

As for the proton injection spectra at all sources, we assume a simple power- 
law spectrum, 

dN 

— oc E- a Q (£ max -E)B(E- E min ) , (5) 



where a, -E max and E min are a spectral index, the maximum acceleration energy 
of protons and the minimum energy of protons at sources respectively. These 
are treated as free parameters. If UHECR sources have a spectral index steeper 
than 2, the total injection energy from a source could exceed the energy budget 
of possible candidates of UHECR sources. Thus, .Emm is required as a cutoff. 
For conserving energetics, a broken power-law spectrum, for example, a = 2 at 
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lower energy than an energy Eq, which is a break energy of the spectral index 
as parameter, was proposed [40,41]. E min in this study should be interpreted 
as a minimum energy if UHECR injection spectrum can be expressed by a 
single spectral index. 

Taking these source-evolution models into account, the total neutrino flux 
from all sources is calculated as 



where / ev (z) and dnj dt(E u , E p , z) are the source-evolution factor and the num- 
ber of cosmogenic neutrinos with energy E v at the Earth produced by protons 
with energy E p injected from a source at redshift z per time at z — 0. The 
normalization factor of the neutrino flux and the value of a are determined 
by fitting the observed UHECR spectra. 



3 Results 

In this section, the results of our calculations of the fluxes of cosmogenic 
neutrinos are described. First of all, we compare our calculated spectrum to a 
spectrum estimated by Ref. [13] to check our simplified treatment of photopion 
production and to investigate the difference of the neutrino spectra predicted 
from 2 different IR/UV background models in Section 3.1. In Section 3.2, the 
normalization factors of the neutrino flux and the spectral indices of UHECR 
injection at sources are determined for every source-evolution model and tran- 
sition scenario so that the calculated spectra of UHE protons best reproduce 
the observed spectra. Then, we investigate the parameter dependence of the 
neutrino fluxes and discuss the capability to judge which transition scenario 
is favored in Section 3.3. The detectability of cosmogenic neutrinos is also 
discussed. 

3.1 Comparison with Another Estimation 

Fig. 2 shows a + spectrum calculated in this study (red solid line) and 
that in Ref. [13] (blue solid line). In addition, our spectrum is divided into 
two spectra: the contributions of the CMB photons (red dashed line) and the 
others (red dot-dashed line). Differences between the two calculations are the 
treatment of photopion production and IR/UV background models. Photopion 
production is treated by a simplified method explained in Section 2 in this 



dN 
dEdt 



(E v ) = J dzf cv (z) J dE p -^(E„,E p ,z) 



(6) 



o 
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Fig. 2. Comparison of our i/„ + v„ spectrum (red solid line) to that calculated in 
Ref.[13] (blue solid line). The same models and parameters are adopted for both cal- 
culations except for the treatment of photopion production and IR/UV background 
models. Our result is divided into contributions by only the CMB (red dashed line) 
and that by the others (red dot-dashed line). The flux difference at ~ 10 17 eV results 
from different IR/UV background models. 

study, while fully simulated by using an event generator SOPHIA [42] in Ref. 
[13]. For IR/UV background photons, we adopt a model constructed by Ref. 
[15], whereas Ref. [13] used another model by Ref. [16]. The other parameters 
are the same ones for comparison: -E max = 10 2L5 eV, E min = 10 18 eV, the proton 
injection spectrum of dN/dE oc E~ 2 - 5 exp (— E/E mSjX ) Q(E — E m[n ), and a 
strong source evolution model denned in Ref. [13]. These fluxes are normalized 
by the same method as in Ref. [13]. 

Our calculation well agrees with the estimation of Ref. [13], but the neutrino 
fluxes at ~ 10 17 eV and above 10 19 ' 5 eV are about a factor of 2 smaller than 
those by Ref. [13]. 
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The discrepancy of the flux at ~ 10 eV is explained by the difference of 
adopted IR/UV background models as follows. At z — 0, the number den- 
sity of photons in the background model of Ref.[15] is about twice as small 
as that of Ref.[16] at ~ 0.2 eV. Since photopion production dominantly oc- 
curs through the delta resonance, s ~ 1.6 GeV 2 , the energy of protons which 
dominantly interact with photons of 0.2 eV is ~ 2 x 10 18 eV and the energy 
of produced neutrinos is ~ 10 17 eV. The mean free path of these protons 
for photopion production is ~ 10 5 Mpc at z = (see Fig. 1), and therefore 
the Universe is transparent against these protons. Thus, the number of inter- 
actions is proportional to the photon number density, and then the IR/UV 
background radiation model by Ref. [16] leads to the neutrino flux twice as 
large as that in the background model used in this study at ~ 10 17 eV . Note 
that these neutrinos are not generated at high-redshift Universe because the 
energy-loss length of Bethe-Heitler pair-creation is much shorter than the in- 
teraction length of photopion production (see Fig.l). Hence, a discussion above 
at z = is sufficient. The shape of the spectrum of cosmogenic neutrinos gives 
us useful information on IR background photons. 

The flux difference above 10 19 5 eV originates from our simplified treatment of 
photopion production. Our treatment predicts a little lower flux than in Ref. 
[13] at the highest energy. 



3.2 Normalization of the Neutrino Flux 

The calculated fluxes of cosmogenic neutrinos are normalized so as to give a 
good fit to observed UHECR spectra in this study. The spectral fit constrains 
the spectral indices of the proton injection spectra and the normalization 
factors of the neutrino fluxes through the normalization of UHECR fluxes at 
the same time. We adopt UHECR spectra observed by the Akeno [43,44] and 
AG AS A [45], and by the HiRes-I and HiRes-II [46] for the normalization. Note 
that the AGASA spectrum is shifted by 10% to lower energies to give a good 
agreement with the Akeno spectrum [47]. 

Fig. 3 represents calculated proton spectra which are best fitted to the observed 
spectra for different source-evolution models listed in the figure. In the upper 
two panels, the calculated spectra are fitted using the chi-square method in 
the range from 10 18 to 10 19 6 eV, which correspond to the proton-dip scenario. 
As a result, spectral indices which best reproduce the Akeno- AGASA results 
are a = 2.7, 2.5, 2.5 and 2.3 for the UNF, QSO, SFR, and SFR+GRBMAC 
source-evolution models, respectively. Fitting to the HiRes spectra leads to 
a = 2.6, 2.5, 2.5, and 2.3, respectively. These indices are independent of E m i n 
and -E max as long as E min < 10 18 eV and E mSuX > 10 20 eV. In the lower two 
panels, we fix the spectral indices to 2.0 and try to fit the calculated spectra 
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Fig. 3. Proton spectra fitted to the Akeno-AGASA data (left two panels) and HiRes 
data (right two panels). Note that a spectrum by the AGASA is shifted by 10% to 
the lower energy to give a good agreement with the Akeno spectrum. £? max and E m \ n 
are set to be 10 22 eV and 10 18 eV, respectively. The upper two panels correspond to 
the proton-dip scenario. The spectral indices of proton sources are 2.7, 2.5, 2.5, and 
2.3 for UNF, QSO, SFR, and SFR+GRBMAC source-evolution models, respectively 
in the upper left panel. In the upper right panel, the spectral indices are 2.6, 2.5, 
2.5, and 2.3, respectively. These are the best fit spectral indices of \ 2 fittings with 
the observational spectra from 10 180 to 10 19 ' 6 eV. In the lower two panels, which 
are in the case of the ankle-transition scenario, the spectral indices are set to be 
2.0. 

to the observed ones in the energy range from 10 19 ' 5 to 10 19 9 eV to represent 
the ankle-transition scenario. In both scenarios, shortfalls in the fluxes at lower 
energies are thought to be compensated by GCRs. Throughout this paper, we 
adopt the normalization based on the Akeno-AGASA spectrum. 

The PAO also reported the energy spectra of UHECRs [48]. However, an 
energy spectrum derived from its ground-based detector covers only down to 
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Fig. 4. Spectra of cosmogenic neutrinos with E m \ n = 10 16 {black) 10 17 {red) and 
10 18 eV {blue) in the proton-dip scenario {left) and the ankle-transition scenario 
{right). The solid lines, dotted lines, and dot-dashed lines show the spectra of z^+z^, 
v e + z7 e , and v e from neutron beta decay, respectively. The cosmological evolution of 
UHECR sources and neutrino oscillation are not taken into account. All fluxes are 
normalized by the Akeno-AGASA spectrum. 

10 18 4 eV. An energy spectrum constructed from a hybrid method covers above 
10 18 eV, but it has a large statistical error due to the small number of detected 
events above 10 19 eV. Thus, we do not use the PAO spectra in this study. 

3.3 Neutrino Fluxes 

We start by investigating the dependence of the fluxes of cosmogenic neutrinos 
on E min . Fig. 4 shows the calculated spectra of + {solid lines), v e + v e 
{dotted lines), v e from neutron beta decay {dot-dashed lines) for E m i n = 10 16 
{black), 10 17 {red), and 10 18 eV {blue). These spectra are calculated based on 
the proton-dip {left) and ankle-transition {right) scenarios. The cosmological 
evolution of UHECR sources and neutrino oscillation are not taken into ac- 
count. The v e + z/ e fluxes are twice as low as the + z/ M fluxes in the energy 
range where neutron beta decay does not contribute to i> e fluxes because pion 
decay produces two muon neutrinos and one electron neutrino. 

An intriguing feature is the spectral peaks of both u^ + u^ and z/ e + z/ e spectra at 
~ 10 16 eV in the left panel. These peaks are generated by interactions between 
protons with ~ 10 17 eV and UV photons. Since there are many protons at 
lower energies due to a steep injection spectrum, UV background photons 
significantly contribute to the neutrino flux. On the other hand, we cannot 
see a spectral peak at ~ 10 16 in the right panel. The spectral peaks are a 
characteristic feature in the proton-dip scenario. 
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We should take neutrino oscillation into account for discussions of the de- 
tectability of cosmogenic neutrinos, as we can only observe the spectrum of 
neutrinos which suffers from neutrino oscillation during their propagation. 
Neutrinos generated from charged pions and successive muon decay have a 
flavor ratio of u e : : v T = 1 : 2 : 0. Neutrino oscillation changes this ratio 
into ~ 1 : 1 : 1 [49]. The flavor ratio approximately holds 1:1:1 even at low 
energies unless anti-neutrinos from neutron f3 decay are dominated. 

Fig. 5 shows the spectra of cosmogenic neutrinos per flavor, Vi + z/j, taking 
neutrino oscillation into account. -Emm and -E max are set to 10 16 and 10 22 eV, 
respectively. Upper limits of neutrino fluxes determined by several experiments 
and the sensitivities of current/future neutrino detectors are also displayed. 
Theoretically estimated fluxes of diffuse neutrinos from AGNs [58] and GRBs 
[59] are also shown. 

As mentioned above, the proton-dip scenario predicts a spectral peak at ~ 10 16 
eV, while the peak does not appear for the ankle-transition scenario. The flux 
difference at 10 16 eV between the two scenarios is about an order of magnitude. 
Therefore, the detection of this peak could be an evidence of the proton-dip 
scenario. The flux of the peak depends on source-evolution models. The SFR 
and SFR+GRBMAC source-evolution models result in three and four times 
higher neutrino fluxes than the UNF source-evolution model, respectively. In 
these strong source-evolution models, the predicted fluxes are comparable with 
the IceCube 3yr sensitivity. 

The neutrino flux at ~ 10 20 eV could be also a clue to judge a better tran- 
sition scenario. The ankle-transition scenario leads to the neutrino flux much 
higher than the proton-dip scenario at the highest energy because of a harder 
injection spectrum. In the SFR source-evolution model, the estimated flux 
reaches the full-ANITA sensitivity. Remember that our estimation of the flux 
of cosmogenic neutrinos above 10 19,5 eV might be underestimated by about a 
factor of 2 as shown in Section 3.1. This works positively for the detection of 
the highest energy neutrinos. The full-ANITA is also expected to detect the 
highest energy neutrinos given the SFR+GRBMAC model is realistic enough, 
and it implies that GRBs are UHECR sources, though we should keep in mind 
that the models of the GRB rate history include large uncertainty [fj. 

The neutrino flux at the highest energy strongly depends on E max . Fig. 6 shows 
the spectra of cosmogenic neutrinos calculated for E max = 10 22 ' , 10 215 , and 
10 21,0 eV in the SFR source-evolution model. All of -E ma x can reproduce the 
observed cosmic ray spectra sufficiently. We can see that the neutrino fluxes 
above 10 19 eV are sensitive to E max . As -E max becomes smaller, the neutrino 

2 A preliminary upper bound of the ANITA approaches the predicted flux of cos- 
mogenic neutrinos with E m3iX = 10 22 eV in the ankle-transition scenario at ~ 10 20 
eV (see \http://www. slac. Stanford, edu/econ f/ C070730/) . 
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Fig. 5. Predicted spectra of cosmogenic neutrinos per flavor (vi+vi) in the proton-dip 
scenario {solid lines) and ankle-transition scenario [dotted lines). These fluxes are 
normalized by using the Akeno-AGASA spectrum. E m { n and -E max are set to 10 16 and 
10 22 eV. The red lines, green lines, and blue lines are neutrino spectra for the UNF, 
SFR, and SFR+GRBMAC source-evolution models, respectively. The fluxes of the 
atmospheric neutrinos (Atm v vertical/horizontal) [50] are represented. As upper 
limits of neutrino fluxes by several experiments, AMANDA-II limits [51,52], limit 
on tau neutrinos by the PAO [53] , ANITA-lite limit [54] are shown. As estimated or 
projected sensitivities, sensitivity of 3 years observation by IceCube [55], sensitivity 
of 5 years observation by the PAO [56], ARIANNA sensitivity [57], and full- ANITA 
sensitivity [54] are also shown. As diffuse neutrino spectra from energetic sources, 
a maximal neutrino flux from active galactic nuclei including neutrino oscillation 
[58] and neutrino spectrum from GRBs calculated in Ref.[59], considering neutrino 
oscillation, with their parameters of £j e t = 1.24 x 10 51 erg, E s h = 10 51 erg, £b = 1> 
£ acc = 100, T = 10 2 - 5 , r = 10 13 - 10 14 ' 5 cm and I = r/Y = 10 10 - 5 cm, which are used 
in Ref.[60], are shown. 

flux at the highest energy is smaller and the flux difference between the two 
scenarios also becomes smaller. In order that the neutrino flux at ~ 10 20 eV 
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can be a clue to distinguish between the two scenarios, sufficiently large E mSLX 
is required. 

Whereas the neutrino fluxes at ~ 10 18 eV are independent of not only the 
transition scenarios (see Fig.5) but also E min and -E max (see Figs.4 and 6), the 
neutrino fluxes are sensitive to source-evolution models. Thus, the neutrino 
flux at ~ 10 18 eV has information on the cosmological evolution of UHECR 
sources. The ARIANNA and PAO will detect cosmogenic neutrinos as long 
as the SFR or SFR+GRBMAC source-evolution scenarios are good ones. We 
can also obtain indirect knowledge on UHECR sources by comparing source- 
evolution models constrained by the neutrino observatories to the cosmological 
evolution models of various astrophysical objects. 

The spectral features of cosmogenic neutrinos described above might be cov- 
ered by the neutrino background from powerful objects like GRBs and AGNs, 
as shown in Fig.5. If the neutrino fluxes from these energetic objects are larger 
than the flux of cosmogenic neutrinos, it would be difficult to test transition 
scenarios by cosmogenic neutrinos. 

A diffuse neutrino flux from GRBs estimated by Ref. [59] is larger than the 
predicted fluxes of cosmogenic neutrinos up to ~ 10 17 eV and hides the spec- 
tral peaks of cosmogenic neutrinos at ~ 10 16 eV which is predicted in the 
proton-dip scenario. However, since neutrinos from GRBs are, in principle, 
distinguishable from cosmogenic neutrinos by time and spatial correlations 
between prompt 7-rays and neutrinos, the spectral peak can be detected. 

For a diffuse neutrino flux from AGNs, the prediction of a maximal contribu- 
tion model by Ref. [58] is shown in Fig. 5. The estimation covers not only the 
peak of the spectrum of cosmogenic neutrinos at ~ 10 16 eV for the proton- 
dip scenario but also the neutrino flux at ~ 10 18 eV. Unfortunately, spatial 
correlation between emitted neutrinos and photons is not expected for dis- 
tant AGNs, because AGNs are not as bright as GRBs in general. However, 
we should notice that this flux is estimated as a maximum. The total flux of 
diffuse neutrinos starts to be restricted by AMANDA-II [52] and will be con- 
strained more precisely by observations in the near future by detectors such 
as IceCube. Whether the spectral features of cosmogenic neutrinos are not 
covered by AGN diffuse neutrinos and are observable will be determined by 
observations. 



4 Discussion 

In Section 3, we found that the spectrum of cosmogenic neutrinos can be an 
indicator to distinguish the two GCR/EGCR transition scenarios though its 
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Fig. 6. The same as Fig.5, but spectra only calculated for the SFR source-evolution 
model in the cases of different E ma _ x , set to 10 22 (red lines), 10 21,5 (green lines), and 
10 21 eV (blues lines). The solid lines are spectra based on the proton-dip scenario, 
and the dotted lines are those on the ankle-transition scenario. 

capability depends on several parameters on UHECR sources: E min , E mauX , and 
cosmological evolution models of UHECR sources. In this section, we discuss 
other uncertainties on the flux of cosmogenic neutrinos and their detectability. 

The difference between the Akeno-AGASA and HiRes spectra brings the un- 
certainty of the neutrino flux through the normalization of the UHECR flux. 
The UHECR flux observed by the HiRes is twice as low as the Akeno-AGASA 
as shown in Fig. 3 as long as the systematic errors of the determination of 
UHECR energies of both experiments are neglected. The flux difference is re- 
flected to the flux of cosmogenic neutrinos directly, and therefore the difference 
of the neutrinos fluxes between the two normalizations is about a factor of 2. 

A calibration of the energy-scale of UHECR observations could decrease the 
uncertainty of the neutrino flux. A dip calibration method in the proton- 
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dip scenario [2,3] leads to a good agreement between the Akeno-AGASA and 
HiRes spectra assuming that the HiRes spectra are shifted to higher energy 
by a factor of 1.2 and the AG AS A spectrum is shifted to lower energy by a 
factor of 0.9. Since the latter shift has already been performed for giving a 
good agreement with the Akeno spectrum, the normalization of the neutrino 
flux based on the shifted HiRes spectra results in the neutrino flux comparable 
with the Akeno-AGASA normalization if the proton-dip scenario is true. 

Since we have not know which transition scenario is favored yet, the uncer- 
tainty on the neutrino flux which originates from the difference of observed 
UHECR fluxes is maximally a factor of 2. 

The number density of IR/UV background photons also causes uncertainty on 
the flux of cosmogenic neutrinos. We already discussed the difference of the 
neutrino fluxes between different IR/UV background models of Refs. [15,16]. 
Since the number density of infrared photons in the model of Ref. [16] is twice 
as large as that of Ref. [15] at low redshift, the former model predicts the 
neutrino flux twice as large as the latter model at ~ 10 17 cV. Here, we also 
compare the best-fit model, which is adopted in this study, with a low-infrared 
model in the same reference. At z — 0, the number density of photons in the 
range of 0. 005-0. 5eV in the best-fit model is about twice as large as that in the 
low-infrared model. This energy range of photons corresponds to the neutrino 
energy of 4x 10 16 -4x 10 18 eV. Thus, a neutrino flux predicted in the low-infrared 
model is twice as small as that in the best-fit model in this energy range. In 
fact, the CMB photons mainly contribute to the total neutrino flux above a 
few xlO 17 eV. The low-infrared model predicts the neutrino flux twice as small 
as the best-fit model at ~ 10 16 -10 17 eV. Adding the result of the discussion 
in Section 3.1 to the discussion above, the uncertainty of the neutrino flux at 
~ 10 16 -10 17 eV is a factor of 4. 

The difference of the spectral index of UHECR injection also slightly changes 
the neutrino flux. The spectral index has been fixed to 2.0 for the ankle- 
transition scenario in Section 3, but somewhat steeper injection spectrum is 
also allowed as mentioned in Section 1. We check the variation of the neu- 
trino flux resulting from the change of the spectral index. Fig. 7 shows the 
comparisons of the neutrino fluxes predicted from two spectral indices in the 
ankle-transition scenario. For reference, a predicted spectrum in the proton- 
dip scenario is also shown. The flux difference between the two transition 
scenarios becomes smaller for a steeper spectrum in the ankle-transition sce- 
nario, which is about a factor of 3 at ~ 10 16 and 10 20 eV. In that case, the 
distinction between the two scenarios becomes little clear. 

The composition of UHECRs could affect the flux of cosmogenic neutrinos. 
The composition is poorly known as mentioned in Section 1. Here, we con- 
sider only a mixed composition model adopted in Ref. [14]. This model is an 
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Fig. 7. Comparison between cosmogenic neutrino spectra with different spectral 
indices in the ankle-transition scenario. The blue solid line and blue dotted line 
are neutrino spectra with a = 2.0 and 2.2, respectively. The red line is neutrino 
spectrum in the proton-dip scenario for reference. These spectra are calculated in 
the SFR source-evolution model with B mffl = 10 22 eV and E m i n = 10 16 eV. Solid 
black lines are the same as those in Fig. 5 

improved version of the ankle-transition scenario. Fig. 9 of Ref.[14] shows that 
nuclei heavier than protons contribute to the neutrino spectrum mainly below 
10 15 eV. Since the neutrinos below 10 15 eV are covered by atmospheric neutri- 
nos as shown in Fig. 5 and 6, protons make a dominant contribution to the flux 
of cosmogenic neutrinos in an observable energy range in the ankle-transition 
scenario. Thus, the capability to distinguish the transition scenarios is not 
dramatically affected. 

Finally, we revisit diffuse neutrino fluxes from AGNs and GRBs and discuss the 
detectability of cosmogenic neutrinos. In Fig. 5, the diffuse neutrino flux from 
AGNs, maximally estimated in Ref. [58], is shown. This estimation assumes 
that GeV-TeV 7-rays from AGNs are of hadronic origin. This scenario requires 
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high-energy hadrons which can generate neutrinos with energies comparable 
to cosmogenic neutrinos. Therefore, a large flux of high-energy neutrinos is 
predicted. On the other hand, the 7-rays can also be reproduced by leptonic 
origin (for review see [61]). This scenario does not always predict UHE neutri- 
nos because UHECR generation is not required. It is still controversial which 
model is the favorable one. Thus, whether the features of cosmogenic neutrino 
spectrum can be observed or not is related to the origin of high-energy 7-rays. 

GRBs are also possible candidates of the sources of UHE neutrinos. In princi- 
ple, high-energy neutrinos from GRBs can be distinguished from cosmogenic 
neutrinos since such neutrinos correlate with the prompt 7-rays of GRBs. 
However, recently, it has been suggested that neutrinos from low-luminosity 
GRBs can contribute to high-energy diffuse neutrino background [62]. Such 
neutrinos might be difficult to distinguish since the prompt emissions from 
such GRBs are difficult to observe. That neutrino flux can be comparable to 
the flux of cosmogenic neutrinos, ~ 10 4 eVm~ 2 sr _1 s _1 at ~ 10 16 eV, but it 
strongly depends on the local rate of the low-luminosity GRBs, which has 
large uncertainties. 

As discussed above, diffuse neutrinos from several powerful objects could pre- 
vent us from observing the spectral features of cosmogenic neutrinos. It is 
determined by future observations whether neutrinos from these objects can 
hide the spectral features. 



5 Conclusion 

We calculated the fluxes of cosmogenic neutrinos for several plausible param- 
eter sets and discussed the possibility that the neutrino flux helps to judge 
which scenario of transition from GCRs to EGCRs is favored. We found that 
the spectrum of cosmogenic neutrinos has a characteristic peak at ~ 10 16 eV in 
the proton-dip scenario as long as extragalactic protons significantly contribute 
to the observed cosmic rays down to 10 17 eV. The predicted flux is comparable 
with the sensitivity of IceCube when the SFR source-evolution model is con- 
sidered. On the other hand, we also found that the neutrino flux at ~ 10 20 eV 
is much larger in the ankle-transition scenario than in the proton-dip scenario 
if the maximum energy of protons generated at sources is sufficiently high. 
If the SFR+GRBMAC source-evolution model is appropriate one, neutrinos 
with ~ 10 20 eV are expected to be detected by Full- ANITA though the flux 
is highly dependent on E max . These spectral features give us clues to judge 
which scenario is favored unless these are covered by the neutrino background 
from potential neutrino sources like AGNs and GRBs. 

We also found that the neutrino flux at ~ 10 18 eV depends on only the cos- 
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mological evolution of UHECR sources. This indicates that the neutrino flux 
at this energy brings us information on the cosmological evolution of UHECR 
sources. The detection of this feature is feasible by ARIANNA and PAO if 
UHECR sources cosmologically evolve like star formation rate. 

As discussed in this paper, the spectrum of cosmogenic neutrinos depends 
on many unknown parameters: E mSuX , E min , source-evolution models, scenarios 
of the GCR/EGCR transition, the shape of the SED of IR/UV background 
photons. That is why cosmogenic neutrinos are the messengers of not only 
the nature of the EGCR sources but also cosmic background radiation. The 
future detection of cosmogenic neutrinos will provide us a lot of information 
on the Universe related to the highest energy phenomena. 
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